P-Vector  Method  for  Determining  Absolute  Velocity  From 
Hydrographic  Data 


ABSTRACT 

Several  major  techniques  ( Stommel-Schott 
method,  Wunsch  method,  and  Bernoulli 
method)  that  have  bent  developed  to  quantita¬ 
tively  estimate  the  geostrophic  velocity  at  the 
reference  level,  ha  ve  the  sa  me  order  of  dynam  i¬ 
cal  sophistication  (geostrophy.  hydrostatic, 
and  density  conservation.)  From  a  technical 
point  of  vieu\  the  Stommel-Schott  method  is 
an  overdetermined  system  (the  number  of  equa¬ 
tions  is  much  larger  than  the  number  of  ixzri- 
ables),  however,  the  Wunsch  method  is  an 
underdelermined  system  (the  number  of  equa¬ 
tions  is  much  smaller  than  the  number  of  vari¬ 
ables).  Based  on  the  same  dynamical  and  ther¬ 
modynamical  framework,  a  simple,  well-posed 
system  (P-vector  method)  is  proposed  in  this 
study.  Consistmit  with  geostrophy,  the  system 
is  assumed  non-dissipaiive.  The  conservation 
of  mass  and  potential  vorticity  leads  to  the  con¬ 
dition  that  the  velocity  vector  is  perpendicular 
to  both  density  ( p)  and  potential  vorticity 
( q  -  fcip/itz)  gradients,  and  that  the  velocity  can 
be  represented  as  Vfx.  y,  z)  =  r(x,  y,  z)P 
(x.  y,  z).  where  P  =  (Np  x  x  V<j(.  The 

unit  vector,  P.  is  computed  from  the  density 
field,  and  the  parameter  r(x,  y,  z)  is  deter¬ 
mined  by  the  thermal-wind  relation.  Further¬ 
more,  an  error  reduction  scheme  is  also  pro¬ 
posed  in  this  study. 


INTRODUCTION 

The  physical  base  for  calculating  geostrophic 
velocity  from  hydrographic  data  is  the  inte¬ 
gral  form  of  the  thermal  wind  relation 
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where  (m,  v).  (u*  v„)  are  the  geostrophic  velocity 
at  any  depth  z  and  at  a  reference  level  z^,  g  is 
the  earth  gravitational  acceleration,  p  is  the 
water  density,  and  p,  is  the  characteristic  value 
of  the  density.  Here,  the  Boussinesq  approxima¬ 
tion  has  been  used.  The  hydrographic  data  only 
determine  the  baroclinie  geostrophic  current, 
i.e.,  the  second  term  on  the  right  hand  side  of 
(1).  The  reference  velocity  (u*  »»)  still  needs  to 
be  determined.  Two  techniques,  the  Stommel- 
Scott  method  (hereafter  referred  as  SS  method) 
and  the  Wunsch  method  (hereafter  referred  as 
CW  method),  have  been  developed  to  determine 
the  reference  velocity  (w,.  vj).  Furthermore,  the 
Bernoulli  method  (Killworth  1985)  has  been  for¬ 


mulated  to  detemune  the  reference  Bernoulli 
function.  B.  and  in  turn  to  compute  the  pressure 
field  (geostrophic  velocity.)  The  SS  method 
(usually  called  the  p  -spiral  method)  is  to  use 
the  conservation  of  mass,  density,  and  vorticity 
to  estimate  reference  velocities  and  diffusivities 
of  heat,  salt,  and  potential  vorticity  (e.g.,  Stom- 
mel  and  Schott  1977;  Schott  and  Stommel  1978; 
Olbers  et  al.  1985).  The  CW'  method  (usually- 
called  the  inverse  method)  is  to  estimate  the 
reference  velocity  by  balancing  the  fluxes  of 
water  mass  and  several  tracers  into  and  out  of 
a  closed  region  (e.g.,  Wunsch  1978:  Wunsch  and 
Grant  1982).  Both  methods  are  based  on  the 
same  order  of  dynamical  sophistication  and  differ 
from  implicit  assumptions  about  the  scales  of 
oceanic  variability  and  different  definitions  of 
the  smooth  field  to  which  the  dynamic  model 
pertains  (Davis  1978).  From  a  technical  point 
of  view,  the  SS  and  Bernoulli  methods  are  over- 
determined  systems  (the  number  of  equations 
is  larger  than  the  number  of  unknowns),  how¬ 
ever,  the  CW  method  is  an  underdetermined 
system  (the  number  of  equations  is  less  than  the 
number  of  unknowns.)  Can  we  design  a  well- 
posed  system  (the  number  of  equations  equals 
to  the  number  of  unknowns)  and  still  keep  the 
same  dynamical  sophistication  as  the  SS,  CW, 
and  Bernoulli  methods?  To  answer  this  ques¬ 
tion.  we  should  first,  investigate  the  conservation 
principles  associated  with  these  methods. 


CONSERVATION  PRINCIPLES 


As  pointed  out  by  Wunsch  and  Grant  ( 19S2), 
in  determining  large-scale  circulation  from 
hydrographic  data  we  can  be  reasonably  confi¬ 
dent  on  the  assumptions  of  geostrophic  bal¬ 
ance,  mass  conservation,  and  no  major  cross* 
isopvcnal  mixing  (except  when  water  masses 
are  in  contact  with  the  atmosphere).  The  poten¬ 
tial  deasity  of  each  fluid  element,  p.  would  be 
conserved,  i.e., 

V-VP=  0  (2) 

where 


V=i 


fix 


m 


Here  V  =  (u,  v.  w).  is  the  three  dimensional  veloc¬ 
ity',  (x.  y.  z )  is  the  coordinate,  and  z  is  positive 
upward.  At  the  undisturbed  ocean  surface, 
z  =  0.  The  potential  density  (potential  tempera¬ 
ture)  is  the  density  (temperature)  a  parcel  of 
water  would  have  if  brought  adiabatically  from 
its  initial  stale  to  the  standard  pressure  of  1000 
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mb.  After  obtaining  the  temperature  and  salinity 
data,  we  compute  the  potential  temperature  (O) 
using  the  Brvden  formula  (Brvden  1973)  and  the 
potential  density  (p)  using  the  equation  of  state 
(the  temperature  be  replaced  by  the  potential 
temperature)  defined  by  the  Joint  Panel  on 
Oceanographic  Tables  and  Standards  (UNESCO 
1981).  Differentiating  (2)  with  respect  to  z.  using 
the  geostrophie  and  hydrostatic  balance,  and 
including  the  latitudinal  variation  of  the  Cori¬ 
olis  parameter,  we  can  obtain  the  conservation 
of  potential  vorticity  equation: 

V-Vq  =  < ),q^/g  (4) 

where /  =  2fl  sin  a,  is  the  Coriolis  parameter,  fl 
=  7.292  x  lO^s*'.  is  the  angular  speed  of  earth 
rotation.  Use  of fdp/dz  may  induce  a  small  but 
systematic  error  into  estimation  of  potential  vor- 
bcity  (Needier  1985).  The  vector  V  is  perpendic¬ 
ular  to  both  Vp  and  Vq,  therefore, 

V  —  Vq  X  Vp  (5) 

Pedlosky  (1986)  pointed  out  that  the  three- 
dimensional  velocity  field  can  be  determined 
from  knowledge  of  the  density  field  unless  the 
potential  vorticity  and  density  surfaces  coincide. 


FIGURE  1  The  absolute  velocity  ant)  the  intersection  of 


TABLE  1  The  vertical  resolution  ol  the  Levitus  climatological  flata  set 
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He  also  credited  Stommel  and  Schott  (1977)  for 
bringing  up  the  fact  that  the  constant  of  integra¬ 
tion  of  the  thermal  wind  equations  is  deter¬ 
mined  if 

Vq  x  Vp  *  0  (6) 

Killworth  ( 1986)  emphasized  the  importance  of 
using  (p.  q)  space  rather  than  (7,  S)  space  for 
determining  absolute  velocities.  Following  these 
ideas  and  taking  a  more  practical  approach,  a 
new  and  simple  method  can  be  constructed  by 
introducing  a  P- Vector  field  to  determine  veloc¬ 
ity  field  by  density  field. 


P  VECTOR  FIELD 


T  Tnder  the  condition  (6),  a  unit  vector  P 

VqXVp 

IVqXVpl 


(7) 


does  exist  and  is  determined  by  the  density  field 
only.  The  P-vector  lies  on  the  intersection  of  the 
constant  potential  density  and  potential  vortic¬ 
ity  surfaces  (Figure  1).  Therefore,  (6)  is  the  con¬ 
dition  (both  necessary  and  sufficient)  for  the 
existence  of  the  P  vector. 

Climatological  P  field  in  the  North 
Atlantic  Ocean  was  constructed  from  the  climato¬ 
logical  mean  potential  density  field  (p)  which 
was  computed  from  the  Levitus  climatological 
mean  temperature  and  salinity  fields  (Levitus 
1982).  The  Levitus  data  has  l°x  1°  horizontal 
resolution,  and  33  vertical  levels  (Table  1).  We 
use  the  finite  difference  to  compute  the  vertical 
and  horizontal  derivatives.  The  horizontal  grid 
is  l°xl°. 


.  2ira  cosi  .  2na 
AX  =  “l«r_'Aj,=  360’ 


and  the  vertical  grid  is  the  difference  between 
the  two  vertical  levels  as  indicated  in  Table  1, 

A zk-z(k)-z(k+  1)_ 

Here,  o  and  a  are  the  latitude  and  earth  radius, 
respectively.  Although  the  unit  vector.  P.  only- 
carries  partial  information  about  the  velocity- 
field,  it  is  rather  surprising  that  the  horizontal 
component,  Pk  =  (P*,  P,),  catches  the  most 
important  features  of  the  North  Atlantic  general 
circulation.  For  example,  at  the  150  m  depth 
P„  (Figure  2)  resembles  both  the  classical  (Fig¬ 
ure  3)  and  recent  (e.g,  Schmitz  and  McCartney 
1993)  views  of  the  North  Atlantic  circulation 
quite  well:  (1)  anticyclonie  subtropical  gyre  in 
the  western  part  of  the  ocean  between  20°-45° 
N,  (2)  recirculation  cell  on  the  western  side 
(west  of  40°  W)  of  the  subtropical  gyre,  (3) 
cyclonic-anticyclonic  dipole  in  the  area  (30°  W- 
10°  W,  30°  .V);  and  (4)  high  latitude  cyclonic  gyre 
(50°-60°  N,  20°  W-50°  W). 
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Comparing  (7)  with  (5).  we  obtain  a 
relationship  between  velocity  and  P  fields: 

V=rfx.  y.  z)P  (8) 

where  r\x,  y.  z)  is  a  proportionality.  Applying 
the  thermal  wind  relation  to  any  two  different 
depths  z,  and  zK,  a  set  of  algebraic  equations 
for  determining  the  parameter  is  obtained: 

=  (9a) 

Au*.  (9b) 

which  are  linear  algebraic  equations  for  r*  and 
r'*'.  Here 

r"1  =  r(x.  y,  2.). 
and 


fp»  J  <>y 

(10a) 

=  p  *Edz- 

(10b) 

The  integrations  in  (10a)  and  (10b)  were  com¬ 
puted  in  this  paper  by  the  trapezoidal  rule, 

k- 1 

f*yf2)d2=5  2 

J  ,m  ~  j  —  7Ti  (11 


If  the  determinant 


FIGURE  2  North  Atlantic  Ocean  P  held  computed  trom  the  Levitus  temperature  and  salinity 
climatology  at  150  m _ _ _ 
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the  algebraic  equations  (9)  and  ( 10)  have  definite 
solutions  for  (m  £  k): 


p;-> 

At’*.  P,'m' 

pm  ptm) 
pm  pw 


(12b) 


As  soon  as  r*'  is  obtained,  the  absolute  velocity 
field  can  be  computed  by  (8). 


SOURCES  OF  ERRORS  OF  THE 
P-VECTOR  METHOD 

For  a  given  level,  z  =  2„  there  are  (N-l)  sets 
(m  =  1, 2, . .  .k-1,  k+ 1, . .  .N)  of  equations 
(9a,b)  for  computing  /*'.  Here,  AT  is  the  number 
of  total  levels  in  the  density  field  All  the  (N-l) 
sets  of  equations  are  compatible  under  the  ther¬ 
mal  wind  constraint  and  should  provide  the 
same  solution.  However,  due  to  errors  in  mea¬ 
surements  (instrumentation  errors)  and  com¬ 
putations  (truncation  errors),  the  parameter  r1* 
may  vary  with  m.  The  error  estimation  and  elim¬ 
ination  of  the  P-vector  method  will  be  discussed 
in  future  papers. 


FIGURE  3  Classical  view  ol  streamline  pattern  in  North 


To  obtain  a  good  solution,  we  should 
first  understand  the  physical  significance  of  the 
P-vector.  Consider  the  P-vector  at  two  different 
levels.  2  =  2..  and  z  -  zt  (Figure  4a),  and 
define  a  to  be  the  angle  between  Pk  and  the 
x-axis  (Figure  4b), 
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P ,  —  cos  <1.  Pi  =sin  q.  (13) 

The  determinant  (12a)  becomes 


where 


p*<  pr' 


p<*‘  pj”1 


sin(a„) 


(14.) 


indicating  the  3  spiral  turning  angle  between 
two  levels,  and  z“m\ 


MAXIMUM  TURNING  ANGLE 
PRINCIPLE 


The  validity  of  the  3  spiral  method  proposed 
by  Stommel  and  Schott  (1977)  depends  on 
the  turning  angle  a„.  If  o„  is  close  to  0  (almost 


FIGURE  4  Vertical  turning  o!  the  P  vector  (a)  P  vector  at 
two  different  levels,  and  (b)  turning  angle  (aOT)  between  z„ 
and  4.  which  is  the  angle  between  the  two  horizontal 
components  P,‘m  and  PJ*. 


z 


V 


no  turning  angle),  the  &  spiral  method  fails  to 
determine  the  absolute  velocity  from  the  hydro- 
graphic  data. 

Different  from  the  3  spiral  method,  the 
P-vector  method  has  (AM)  choices  (degree  of 
freedom)  for  picking  a  level  (zm)  for  the  compu¬ 
tation.  Theoretically,  all  the  (AM)  sets  of  equa¬ 
tions  (9)  are  equivalent.  Therefore,  we  have  free¬ 
dom  to  ignore  those  levels  leading  to  very  small 
turning  angles.  To  assure  the  best  solution,  we 
adopt  a  maximum  turning  angle  principle  to  get 
the  level 

jsin(ttta,t)i  =  ^^isin(tOj  (15) 

for  computing  r1*-  in  (9)  and  in  turn,  determining 
the  absolute  velocity  V^ix,  y,  a,)  at  the  level  z,. 
The  horizontal  velocity  field  F,<0;  at  the  500  m 
depth  (Figure  5a)  displays  a  circulation  pattern 
in  the  main  therm oc line  very  similar  to  that 
obtained  by  Olbers  et  al.  (1985)  using  SS  method 
(Figure  6).  The  circulation  pattern  is  dominated 
by  a  broad  Gulf  Stream  leaving  the  North  Ameri¬ 
can  coast  at  about  35°  A',  turning  tow  ards  the 
east,  at  3 8®  N,  60°  W'  to  a  branching  point  at  about 
40°  A",  40°  W,  splitting  into  two  branches.  The 
northern  branch  moves  water  northeastward 
until  reaching  50°  AT  and  deflects  towards  north 
at  50°  AT,  15-30°  W.  The  width  of  this  branch 
increases  northeastward  from  700  km  at  40°  N, 
50°  IV  to  nearly  1500  km  in  the  northeast  Atlantic 
at  50-60°  N.  15-30°  IV. 

The  southern  branch  moves  water 
southeastward  between  30°  Ar  and  35°  A'  toward 
Gilbratar.  and  feeds  into  the  subtropical  gyre. 
Our  solution  also  shows  a  tight  recirculation 
cell  on  the  western  side  of  the  subtropical  gyre 
which  was  depicted  in  earlier  research  (e.g., 
Stommel  et  al.  1978;  Reid  1978:  Olbers  et  al. 
1985).  Furthermore,  a  northern  cyclonic  gyre 
exists,  but  is  much  weaker  than  the  subtropi¬ 
cal  gyre. 

Figure  5b  shows  the  vertical  velocity 
at  500  m  depth,  computed  by 

w~r(x,  y,  z)P,(.x,  y,  z). 

In  the  vast  areas  of  the  North  Atlantic  Ocean, 
the  vertical  velocity  is  quite  weak  (less  than  1 
m/day)  and  noisy  except  in  a  very  small  area 
south  of  Iceland,  where  a  very  strong  downward 
velocity  (>  4  m/day)  is  found.  This  coincides 
with  the  strong  downward  branch  of  the  global 
conveyor  belt  (Broecker  1991).  We  also  see  that 
in  the  subtropical  gyre,  the  water  is  usually  down- 
welling.  The  P-vector  method  seems  to  have  a 
capability-  to  diagnose  the  vertical  velocity- 
signal. 
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ERROR  REDUCTION  SCHEME 

W'e  refer  the  velocity  field  F5:  directly  com¬ 
puted  from  the  P-veetor  to  the  Oth-order 
solution.  We  might  ask  if  the  Oth-order  solution 
realistic?  We  need  to  check  the  conservation 
laws  (2)  and  (4).  The  computed  velocity  field 
contains  computational  noise  due  to  the  trunca¬ 
tion  errors  and  data  noise  due  to  the  instrumen¬ 
tation  errors.  We  use  the  iteration  method  to 
eliminate  the  noise.  Let  L  denote  the  iteration 
step.  L  =  0  refers  the  velocity  field  V"n.  Is  V°‘  a 
good  solution?  We  need  to  use  the  local  and 
global  conservation  indices  to  verify. 


Local  and  Global  Conservation  Indices 

The  errors  may  lead  to  a  non-conserva¬ 
tion  of  mass  and  potential  vorticity,  as  mea¬ 
sured  by 


J*®-| cos  ««•>)  = 
K9L)m\ cos  a,,X))  = 


|V^Vp| 

I^IIVpi 

I^IN 


(16a) 

(16b) 


where  a  J,0  is  the  angle  between  VL-'  and  Vp  q 
f 1  the  angle  between  F !t)  and  Vq,  R^Xx,  y,  z) 
and  R^Xx,  y,  z )  the  scaled  residues  of  the  den¬ 
sity  and  potential  vorticity  conservation  equa¬ 
tions  at  each  point  for  the  L-th  iteration.  Theo¬ 
retically.  F  should  be  perpendicular  to  both  Vp 
and  Vq,  i.e.. 


Rv(x,  y.  z)  =  0.  R,(x.  y.z)=  0 

for  all  locations.  We  define  the  local  conserva¬ 
tion  index  (LCI)  as 


R.  LXx,  y,  z)  =  R,"~(x,  y.  z)+R;'\x,  y.  z). 

C17a) 

and  the  global  conservation  index  fGCI)  as  the 
average  of  LCI  over  the  whole  domain 

^  (17b) 

where  .V  is  the  total  number  of  the  data  points. 

Local  Reference  Level 

The  values  of  LCI  represent  the  accu¬ 
racy  of  the  solution  at  each  point  (x.y.z).  How¬ 
ever,  the  value  of  GC1  denotes  the  overall  accu¬ 
racy  of  the  solutions.  The  smaller  the  values  of 
these  indices,  the  more  accurate  the  velocity 
field  would  be.  The  two  parameters  LCI  and 
GCI  have  different  usages.  LCI  is  used  for  pick¬ 
ing  the  local  reference  level  and  GCI  for  neces¬ 
sariness  of  the  next  iteration.  Starting  from  the 
velocity  field  F°  obtained  from  the  P-vector 
method,  we  check  every  water  column  (x„  y,), 
compare  LCI  vertically,  and  choose  the  level 
with  the  minimum  LCI  value 

R*(x,  y.  z^)  =  ™nR,(x.  y.  a)  (18) 

as  the  reference  level  for  that  column.  After 
Zrrfyc,  y)  is  determined,  we  can  use  (1)  to  re¬ 
compute  the  velocity  field.  F which  is  called 
the  first  iterative  velocity. 

Iteration  Method 

Is  V1'  better  than  V"0'?  We  can  use  GCI 
to  compare  V,1>  and  F'0'.  The  condition 

GCP^GCP1  (19a) 

indicates  that  F0)  is  better  than  tA1).  We  do  not 
need  to  do  any  iteration  and  choose  Ft0:  as  the 
optimal  solution.  The  condition 

GCP><GCP»  (19b) 

denotes  that  the  V1,  is  better  than  l'10’.  We  need 
to  pursue  a  new  iteration,  i.e.,  computing  LCP\ 
in  turn  determining  z^1'  for  F111  and  calculating 
F17'  by  using  (1)  and  Comparison  between 
GCP'  and  GCP 11  leads  to  a  choice  if  a  further 
iteration  is  needed  (Figure  7).  Such  an  iteration 
process  repeats  until  no  further  improvement 
will  be  made. 

ABSOLUTE  VELOCITY  AT 
DIFFERENT  DEPTHS 

The  absolute  velocity  vectors  at  different  lev¬ 
els  are  shown  in  Figure  8.  The  flow  pattern 
at  100  m  level  (Figure  8a)  is  quite  similar  as  that 
at  500  m  level:  the  Gulf  Stream,  an  anticyclonic 
subtropical  gyre,  a  weak  northern  cyclonic  gyre. 
Comparison  between  Figure  5a  and  Figure  8a 
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leads  to  the  fact  that  the  subtropical  gyre  shrinks 
m  its  north-south  extension  with  increasing 
depth,  which  agrees  with  Olbers  et  al.'s  (1985) 
results. 

The  flow  pattern  at  deep  ocean  (2500 
m.  3500  m.  and  4500  m  levels)  is  shown  in  Figure 
Sb-d.  The  dominant  features  at  2500  m  level  (Fig¬ 
ure  8b)  are  the  western  boundary  currents. 

After  passing  the  Gibbs  fracture  zone  at  about 
50°  A’,  40°  W  the  current  splits  into  two 
branches:  a  weak  northern  branch  and  strong 
southern  branch.  The  southern  branch  flows 
southwestward.  turns  westward  at  40°  A’  50°  W, 
and  then  follows  around  the  east  coast  of  North 
America  As  this  branch  enters  the  Sagasso  Sea 
Trench,  a  noticeable  cyclonic  eddy  appears  at 
4500-m  level  (Figure  8d)  to  conserve  the  poten¬ 
tial  vorticitv  as  the  ocean  depth  increase.  The 
maximum  swirl  speed  is  2  cms  :.  After  passing 
30°  Ar.  the  mid-depth  flow  at  the  2500  m  and  3500 
m  (Figure  8c)  levels  recirculates  northeastward 
at  30''  .V.  70"  W’  with  a  typical  speed  2  cmsr'. 
Different  from  SS  method,  our  solution  does  not 
showr  a  strong  southeastward  flow  near  Baha¬ 
mas  carrying  water  into  the  South  Atlantic 
Ocean. 

A  cvclonie-anticyclonic  pair  is  found  in 
the  deep  water  at  the  eastern  (20c-40c  W)  tropi¬ 
cal  part  (10°-30°  A)  of  the  North  Atlantic  Ocean. 
The  cyclonic  eddy  is  in  the  south,  and  the  anticy- 
clonic  eddy  is  in  the  north.  The  maximum  swirl 
speed  is  around  2  cm/s.  This  feature  is  indeed  . 
strikingly  similar  to  the  map  (Delant  1941)  of 
the  absolute  flow  at  2000  m  (reproduced  by 
Reid  1981).  Furthermore,  our  computation 
shows  that  this  dipole  structure  becomes  more 
evident  as  the  depth  increases. 

The  absolute  velocity  computed  by 
Olbers  et  al.  (1985)  has  less  noise  at  the  deep 
level  low  latitudes  than  the  P-Vector  method. 
This  is  partly  because  the  turbulence  fluxes  and 
associated  mixing  are  neglected  in  our  computa¬ 
tion. 

ABSOLUTE  VELOCITY  AT 
DIFFERENT  VERTICAL  CROSS 
SECTIONS 

To  compare  the  results  from  this  method  with 
the  geostrophic  calculation  based  on  the 
results  of  the  CW  method  (Wunsch  1978;  Wunsch 
and  Grant  1982)  and  the  SS  method  (Olbers  et 
al.  1985),  we  display  vertical  distributions  of 
absolute  velocity  in  latitudinal  sections:  57°  (V. 
30°  W  (Figure  9)  and  zonal  sections:  24°  A’.  53°  A 
and  59°  A  (Figure  10).  These  figures  present 
the  velocity  component  normal  to  these  sec¬ 
tions. 

We  see  some  main  features  in  the  west¬ 
ern  (Figure  9a)  and  eastern  parts  (Figure  9b) 
of  the  subtropical  gyre:  (a)  a  strong  upper  ocean 


FIGURE  7  The  iteration  method  tor  error  reduction 
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(above  1000  m  depth)  eastw  ard  flow  (the  Gulf 
Stream)  between  30-40°  K  with  a  maximum 
speed  greater  than  20  cm/s,  (b)  a  weak  west¬ 
ward  flow  below  the  Gulf  Stream  (western 
boundary  currents),  (c)  a  near  surface  west¬ 
ward  flow  south  of  30°  A  (southern  branch  of 
the  subtropical  gyre)  with  a  maximum  speed 
of  6  cm/s.  and  (d)  banded  structure  in  the  deep 
layers  (z  <  -1500  m),  where  the  current  direc¬ 
tion  alternates  on  a  horizontal  scale  around  2000 
km,  which  is  wider  than  the  Wunsch  and  Grant 
(1982)  results. 

We  also  see  some  interesting  features 
of  the  meridional  flow  at  different  latitudes 
(24°  A,  36°  A,  53°  A,  59°  A),  which  is  summarized 
as  follows,  (a)  In  the  upper  ocean  (above  1000 
m),  at  24°  A  (Figure  10a)  a  narrow  northward 
flow  appears  near  the  east  coast  of  North 
America  (the  Gulf  Stream)  from  65°  W-75°  W, 


MTS  Journal  •  VoL  29,  No.  2  •  9 


LATITUDE 


FIGURE  8  ABsotuie  velocity  vectors  at  ditlerent  levels:  (a)  100  m.  (bi  25Q0  m.  (c)  3500  m.  anc  (fl)  4500  m 
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with  the  maximum  speed  around  5  cm/s  and  towards  east  as  the  latitude  increases  (Figure 
the  width  near  1000  km.  and  a  weaker  returning  10c  and  d):  the  maximum  current  speed  appears 
southward  flow  occupies  the  rest  of  the  region  between  longitudes  20°-30°  W.  (b)  In  the  deeper 

(65°  W-20°  WO.  At  the  latitude  of  36°  N  (Figure  layer  (below  1000  m  depth),  a  banded  structure 

10b),  a  narrow  southward  flow  appears  revealing  the  alternating  northward  and  south- 

between  60°-70°  N  and  is  sandwiched  between  w’ard  flow  with  smaller  current  speed.  The 
two  northward  flows.  This  indicates  the  exis-  width  of  these  alternating  bands  is  around  1000- 
tence  of  a  recirculation  cell  at  that  latitude  (36°  2000  km,  which  is  also  wider  than  the  results 

N).  The  core  of  the  northw  ard  flow  shifts  of  Wunsch  and  Grant  (1982). 
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FIGURE  10  Continued 
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Comparing  Figures  9  and  10  with  Fig¬ 
ure  18a-d  in  Others  et  al.  (1985)  and  with  the 
results  in  Wunsch  and  Grant  (1982),  we  found 
that  our  results  are  closer  to  the  SS  method. 

CONCLUSIONS 

The  P- Vector  method  provides  a  very  simple 
scheme  to  compute  the  absolute  velocity 
from  hydrographic  data.  The  resulting  circula¬ 
tion  pattern  resembles  the  classical  view  of  the 
North  Atlantic  circulation  (Wust  1935:  Defant 
1941).  especially  in  the  upper  ocean.  The  Physi¬ 
cal  base  of  this  method  is  the  conservation  of 
potential  density  and  potential  vorticity.  There¬ 
fore,  if  we  use  a  synoptic  hydrographic  data  set 
to  determine  the  absolute  velocity,  we  will  get 
better  results  since  a  synoptic  snapshot  of  the 
ocean  fields  obey  the  conservation  laws. 

Similar  results  were  obtained  using  the 
P-Vector  method,  which  does  not  incorporate 
turbulent  mixing  and  fluxes,  and  the  SS  method, 
which  does  use  turbulent  mixing  and  fluxes 
(Olbers  et  aL  1985).  Both  use  the  same  Le virus 
data  set.  This  indicates  that  the  main  features 
of  the  North  Atlantic  Ocean  circulation  are  deter¬ 
mined  by  the  density  field  rather  than  by  turbu¬ 
lent  fluxes  and  mixing.  This  supports  Pedlosky's 
idea  that  the  density  field  itself  can  determine 
the  main  features  of  the  ocean  circulation. 

The  dependence  of  geostrophic  cur¬ 
rents  on  only  the  density  field  will  lead  to  a 
better  initiali2ation  scheme  for  ocean  numerical 
models.  Ocean  numerical  models  are  generally 
initialized  with  either  observed  or  climatological 
temperature  and  salinity,  as  well  as  zero  veloci¬ 
ties.  This  leads  to  an  imbalance  between  the 
velocity  and  density  fields  during  the  spin-up 
period.  If  we  use  the  geostrophic  velocity,  com¬ 
puted  by  P-Vector  method  (from  the  density- 
field  only)  instead  of  zero  velocity,  as  an  initial 
condition  for  the  numerical  models,  the  velocity 
and  density  fields  are  balanced.  This  might 
reduce  the  spin-up  time. 
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